Spin waves in quasi-equilibrium spin systems. 
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Using the Landau Fermi liquid theory we have discovered a new regime for the propagation of 
spin waves in a quasi-equilibrium spin systems. We have determined the dispersion relation for the 
transverse spin waves and found that one of the modes is gap less. The gap less mode corresponds to 
the precessional mode of the magnetization in a paramagnetic system in the absence of an external 
magnetic field. One of the other modes is gapped which is associated with the precession of the 
spin current around the internal field. The gapless mode has a quadratic dispersion leading to some 
interesting thermodynamic properties including a J 13 / 2 contribution to the specific heat. We also 
show that these modes make significant contributions to the dynamic structure function. 

PACS numbers: 



There has been much renewed interest in studying non- 
equilibrium (NEQ) spin polarized paramagnetic Fermi 
liquids and gases. The NEQ spin system can be de- 
fined as a polarized spin system which can be ob- 
tained by maintaining an unequal population of differ- 
ent spin species to create an off-balanced Fermi level in 
the absence of an external field. This is both an old 
subject 0,0], and a rather new one[E 0, IE IE 0- The 
NEQ spin systems have been created in a wide variety of 
materials using many novel techniques. An example of 
the most recent one would be the creation of steady state 
spin polarization in a system of cold and ultra-cold Fermi 
gases of alkali metals. In these system one uses optical 
trapping and atomic manipulation to get unequal popu- 
lations of the fermions in two different states [EGl ■ The 
NEQ spin systems have also been created in many mate- 
rials which comprise the study of one of the most promis- 
ing areas of physics, namely spintronicsjE IE IE EJ- I n 
this area of study a couple of methods are in practice to 
achieve the NEQ state of these spin systems. 

Widely used methods for achieving NEQ states are op- 
tical pumping and polarized spin injection. In the opti- 
cal pumping method 6], electron orbital momentum is 
directly oriented by a strongly circularly polarized laser 
light and through the spin-orbit interaction the elec- 
tron spin is polarized. In the polarized spin injectionjfj 
method one connects a magnetic electrode to a normal 
metal sample and drives a polarized current using ex- 
ternal biasing which creates unequal population of spins 
in the normal metal. In such a ferromagnet - metal bi- 
layer structure, it has been suggested that if the direction 
of the magnetization in the ferromagnet is made to pre- 
cess (by applying an external magnetic field on the fer- 
romagnet) we can inject polarized spin into the normal 
metaipjj]. It has also been suggested in the literature 
that we can use the interaction of the electron spin with 
the strain-induced potential to create the NEQ system 
in a nano-mechanical system The NEQ system of po- 
larized liquid helium is one example of an old but still 
interesting subiect|ll|. The rapid melting method sug- 
gested by Castaing and Nozieres [L2l fl3| . and the optical 
pumpingp], |E El were widely used to create the NEQ 



systems of liquid helium. 

In this Letter we are interested in studying the ther- 
modynamics and the spin dynamics of the NEQ systems 
in the limit of small spin polarization. We know that 
for a paramagnetic spin system of free fermions in the 
presence of an external magnetic field, the magnetiza- 
tion, m, precesses around the applied magnetic field at 
the Larmor frequency. For a small transverse perturba- 

— * 4- 

tion 5 B = B x x + B y y, and defining SB = B x ±iB y and 
tyv^ = m x ± im v , the evolution of m for non interacting 
Fermi system can be written as|l4j . 

dm* , . , 

— — = 2im ± jB (1 
dt ' w 

where m = moe~ ia, ° * and uj^ = ±2jB is the Larmor 
frequency and B is the external magnetic field. The 
fact that the magnetization precesses around the exter- 
nal field with the Larmor frequency also holds true for 
a system of interacting fermions, with spin conserving 
interactions. Now, for the NEQ systems, since the exter- 
nal field is zero, the precessional mode should be gapless. 
To understand the details of the dispersion of the gap- 
less mode we need to study the thermodynamics and the 
dynamics of a Fermi liquid away from equilibrium. 

The general theory of Fermi liquids and gases far 
from equilibrium is beyond our current understanding. 
We will develop here the theory for weakly spin polar- 
ized paramagnetic Fermi liquids with B = 0. In this 
limit we can argue that the Fermi liquid will be in a 
quasi-equilibrium(QEQ) state. It is well known that 
the polarization dependent chemical potential has the 
form fi(m) — fi(0) + Q(ro 2 )|l5||. hence for weak polar- 
ization, i.e. for - <C 1 where n is the particle density, 
iu(m) f» m(0). This suggests that for small m the system 
will be in a QEQ phase and we can study the dynamics 
of the spin system around this phase. We know that the 
variation of the energy per-unit volume due to a varia- 
tion in the distribution function from the ground state is 
written as[TllTllT7|: 

e = £ lJ n p<y + fpp' 5n pa 5n p i a . + .... (2) 

P<7 p(J,p'(j' 



2 



So the quasi-particle energy can be written as: 
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where for a spin conserving interaction with ^ < 1, 

£<j<j' ^ rs i fa f 
J ppi — J ppi ~r J ppi u ' o 
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Using Eq.@ and expanding in Legendre polyno- 
mials we can derive that the chemical potential of the up 
and the down spins will have the forms /i = e F ^ + jf^S- + 

0(m 2 ) and n = e a Fl -§^ + 0(m 2 ), where N(0) = ^f- 
is the m = density of states of the quasi-particles at the 



Fermi level. So in leading order in m, e F + - 
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For the EQ system, fi = e^T — B + + (9(m 2 ) and 



N(0) 1 

= 2B 



0(m 2 ). 
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;ain in leading order 
Contrary to the 
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in m, Ef-\ — £fi 
EQ system, we see that the change in the Fermi energy 
of the up spin and the down spin fermions of the QEQ 
system depends explicitly on the Fermi liquid parameter 
Fq . For the QEQ and EQ systems with a majority of 
up spins, we have shown, in Fig. la, the change in the 
Fermi energy of the up spin (ep^) and the down spin 
(e'pi) fermions, and in Fig. lb, the energy spectrum of 
these spin species. Here B e ff is the effective field. For 

the QEQ system B e jf — ~^q° and for the EQ system 
Beff — wgy = iqrp^- F° r a QEQ system, Fig. la cor- 
responds to a steady state with constant number den- 
sity which can be obtained, for example, either by taking 
out some down spins and putting in equal number of up 
spins or by flipping the down spin over. To maintain such 
steady states we need to continuously supply energy to 
the system (e.g., continuously shining light in the system 
in the optical pumping method) against the loss in the 
polarization due to spin flip scattering mechanism. In 
what follows we want to fix the number density of the 
Fermions to avoid any complication that would be intro- 
duced due to the coupling between the charge density 
fluctuations and the spin density fluctuations. 

To study the collective excitations of the QEQ sys- 
tem we investigate the oscillations of the transverse com- 
ponent of the magnetization (Sm p ). The linearized ki- 
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netic equations for the evolution of Sr 
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_> -> _> dn° p -> 
v p ■ V(f5m p - -t^qS h p ) 
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where Sh p = —SB 



J2 P ' fpp'^ m p' i s the fluctuation in 
the effective magnetic field and 7[ TO p] is t ne collision in- 
tegral. The equilibrium field h® = —B + J2 P < fp P ' m p' ■ 
For the QEQ system B — . We can use the re- 
laxation time approximation for the collision integral 



FIG. 1: For the QEQ and EQ systems with a majority of up 
spins: (a) the change in the Fermi energy, and (b) the energy 
spectrum of the up and the down spin fermions. -Be// = 
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for the QEQ system and B e ff = jjmj = 



1+Fg 



for the 



pi — 

_( "»*-■"» ) where to is the spin diffusion 
lifetime where in a Fermi liquid it varies as T~ 2 with the 
temperature. Here we investigate the collective modes of 
the QEQ system at T — > OK for which we can set 4 TO p] 
to be zero. We will discuss briefly the effect of the finite 
temperature on the modes later in this paper but the de- 
tails of the calculation will be shown elsewhere. We have 
also assumed that nuclear relaxation time, 7\, is suffi- 
ciently long so that we can maintain the QEQ state for 
time much longer than any diffusion time. The effect of 
the finite T\ is being studied currently. 

Defining m p = -|jr^y , and Srn^ = -j^rVp HEEll 

and using v^(r,t) = j d 3 qdujvt-{q, uj)e 1 ^ q ' r ~ wt ) for the 
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To derive dispersion relations we set SB 



0. We 



substitute v p by its expansion in Legendre polynomials: 

Vp = J2i v^PiiP ' o) an d project out the I = 0, I = 1, 
I = 2 and I = 3 components of Eq.©. Because of the 
cos 9 term in ■ ~q one sees that every moment couples 
to one higher moment and one lower moment. Solving 

the equation for I — 3 moment we find that 

which can be negligible for low q. We also find that 
the exclusion of the I = 2 moment also does not change 
the important aspect of the dispersion. To see the main 
features of the modes we will first work with the pro- 
jection of the I — 0, and I = 1 moments and truncate 
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the distortions of the Fermi surface, 
parameters, f", for I > 2. 

For the I = moment, we get, 



and the Landau 
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^qvpvf = 



(7) 



which gives the net magnetization conservation law. 

For the I = 1 moment, which gives the motion of the 
magnetization current, we get, 



(w + ojf)vf - qv F (l + Fg)i 
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FIG. 2: Normalized ui for the collective modes and the p-h 
continuum as a function of normalized q for the QEQ spin 
system. 
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FIG. 3: Normalized u for the collective modes and the p- 
h continuum as a function of normalized q in the EQ spin 
system. 



where tuf = ±^j^(Fq is related to the internal 

field due to the interacting fermions. 
Solving the Eos.l|7l8|>. we get, 



ch 2 



(16) 



lu(lu + ujf) - c 2 s q 2 = 



(9) 



where <? s = + Ffi)(l + ^f)v 2 F . We solve Eq.Q for u. 
The solutions are: 



2„2 



U\ (q) = + 



(10) 



(11) 



If we take Fj 1 and = for I > 3 and include the 
projection of I = 2, Eq.@ will be modified to, 

( W ( W + u>f) c 2 s q 2 )(cu + w ±) - % 2 sl q 2 = (12) 

5 

where c 2 sl = ±(1 + ^)(1 + ^f-)v 2 F . From Eq.JTSJ we can 
show that the contribution of order q 2 to the dispersion 
relations will be only to u>^(q), which now reads as, 



uf{q) =ujf + 



c 2 q 2 AN{Q)vl, 3 . 2 

T -J- T —^( h l)q 

30m y F? H 



(13) 



For the EQ system, Eq.© will have the form 
(wTwo-V9± 2fgm)i^ - {v* v • ~q N(0) ± 2m) 
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which modifies the Eas. H10|) and i|ll|) as follows: 
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where 



N(0) 
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B. We mention here, for the purpose of 



later use, that the particle- hole (p-h) continuum for the 
QEQ and the EQ systems are given by uj^ h = =F2m/g + 

~q.~v and ui^ h — ±cj T 2m/ a + .~v , respectively. 

We evaluated the dispersion relations of the collective 
modes and the p-h continuum for 10% polarized QEQ 
and EQ spin systems using the Fermi liquid parameters 
suitable for a polarized liquid helium at zero pressure [l^. 
For the EQ system ^ = 0.04 corresponds to 10% polar- 
ization. The results have been presented in Figs. (2,3). 

In Fig. 2, we see that the QEQ system has one gapless 
and one gapped collective mode. In principle, we will 
have / gapped modes if we include up to I th term of the 
expansion of f££! and in the Landau kinetic equation. 
We have shown one of those gapped modes in Fig. 2. The 
gapless mode has a dispersion relation similar to that of 
the Goldstone mode of a ferromagnetic system^)]. The 
gapless mode is related to the precessional mode of the 
magnetization in the absence of an external field. This 
fact can be understood by setting q — in Eq.J7J), which 
gives uiv^ — 0, and for ^ we get ui = and hence 
a gapless mode. The gapped mode originates from the 
oscillation of the spin current, j a oc v\. From Eq.JSJ with 
q = we have (u> + uf)^ = 0, and for i/^ ^ we 
find oj = ojf. Thus the gapped mode is related to the 
precession of the spin current around the internal field 
due to the interacting fermions. We also see that these 
collective modes are well separated from the p-h contin- 
uum at low q and hence will be propagating modes. If 
we did not include /f in the kinetic equation, the cur- 
rent mode would lie inside the p-h continuum, and hence 
would be Landau damped. This shows the importance 
of including in the kinetic equation. In Fig. 3, we see 
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that both the collective modes of the EQ system are also 
propagating. But contrary to the QEQ system, none of 
the modes are gapless. The uniform mode in the EQ sys- 
tem corresponds to the precession of the magnetization 
around the external field at the Larmor frequency ±loq. 
The current mode corresponds to the precession of the 
spin current and is gapped with iujQ + uff . 

The presence of the /f in the kinetic equation has 
an interesting consequence in the dynamical structure 
function, S(q,u>). From Eq.©, we find an expression 
for the dynamical spin susceptibility using x ± (q,w) = 

=^^ 01 03 and used ii; t0 obtain 

S ± (q,u>) = A/m^faw)]- Jm^fa-a;)]) (17) 

Z7T 

to study the frequency dependence of S^(q,ui) and the 
f-sum rule. For experimental purpose S^(q,uj) is de- 
fined in relation to the scattering cross section as Sjjf- — 
aoS^lq, u)\2(ty where <7o contains all the elements related 
to the probe particle interaction and S ± (q, u>) character- 
izes the many-body system. The f-sum rule states that 
/ ojS ± {q,oj)duj = (1 + £l)3EL [23| , where n is the to- 
tal density, and m* is the effective mass. For the QEQ 
system we see that the gapless and the gapped collective 
modes exhaust the sum rule if we include the /f in the 
kinetic equation. In this case the p-h contribution to the 
sum will be proportional to q A . But if we did not include 
fi the p-h continuum would also have a contribution pro- 
portional to q 2 . The gapped mode, which comes out of 
the p-h continuum due to the inclusion of ff in the ki- 
netic equation, takes out most of the spectral weight of 
the p-h continuum contribution to the sum rule. We ob- 
serve a similar behavior for the EQ phase, as was found 
earlier in the ferromagnetic system[Ljj. 



In addition to the dynamical consequences of the gap- 
less mode there may be interesting thermodynamic prop- 
erties that can be observed in the QEQ system. Perhaps 
the most dramatic would be a T 3 / 2 contribution to the 
specific heat. Whereas, for the EQ system, because all of 
the transverse collective modes are gapped, the specific 
heat contribution would decreases exponentially with the 
decrease in temperature. 

As mentioned earlier, to study the effect of the tem- 
perature on the collective modes, we solved the Landau 
kinetic equation using the relaxation-time approximation 
for the collision integral. We find that the real part of 
the dispersion relations remains dominant at low temper- 
ature. The dispersion relation shows a transition of the 
spin system from the collision less region to a hydro dy- 
namic region as has been discussed in the literature |16(. 

In summary using the Landau Fermi liquid theory we 
explored the thermodynamics and the dynamics of a 
paramagnetic spin system which is out of equilibrium. 
We report for the first time the prediction of a gapless 
collective mode in such systems. We discuss the differ- 
ence in the nature of the collective modes in the QEQ 
and the EQ systems. We point out the importance of 
including the /f Landau parameter in the kinetic equa- 
tion. We believe that some of the qualitative features of 
the collective modes, in particular the gapless mode, will 
also be seen in systems far from equilibrium. We want to 
emphasize that this study is applicable in a wide variety 
of interesting materials. 
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